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In nucleus-nucleus collisions, high-pr electron spectra depend on the medium mod- 
ified fragmentation of their massive quark parents, thus giving novel access to the 
predicted mass hierarchy of parton energy loss. Here we calculate these spectra in a 
model, which supplements the perturbative QCD factorization formalism with parton 
energy loss. In general, we find - within large errors - rough agreement between theory 
and data on the single inclusive electron spectrum in pp, its nuclear modification fac- 
tor R'xai ^iid its azimuthal anisotropy t"!- However, the nuclear modification factor 
depends on the relative contribution of charm and bottom production, which we find 
to be affected by large perturbative uncertainties. In order for electron measurements 
to provide a significantly more stringent test of the expected mass hierarchy, one must 
then disentangle the b- and c-decay contributions, for instance by reconstructing the 
displaced decay vertices. 



1. Introduction 



The semi-leptonic decays of charmed and beauty mesons dominate the high-p^^ 
electron spectrum in ^/sNN = 200 GeV hadronic collisions at the Brookhaven 
Relativistic Heavy Ion Collider (RHIC) up to pt ~ 20 GeV, where the Drell- 
Yan contribution starts to become significant. At the 30 times higher center 
of mass energies of the CERN Large Hadron Collider (LHC), heavy quark de- 
cays are expected to dominate the electron spectrum up to pt ~ 30 - 35 GeV, 
where I^-decay contributions take over. In relativistic nucleus-nucleus colli- 
sions, such electron spectra are of great interest, since they are expected to 
give access to the medium modified fragmentation of their heavy quark par- 
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ents [1-4]. Models of strongly medium-enhanced fragmentation of light quarks 
and gluons [5-9] are favored by RHIC data, which show a strong suppression 
of high-pT hadron spectra and hadron-triggered correlation functions [10-13]. 
Establishing the dependence of this suppression on the mass and color charge 
of the parent parton provides a novel opportunity to further test the micro- 
scopic picture conjectured to underlie medium-induced high-pr hadron sup- 
pression [14]. This may also help to further constrain information about the 
density and collective dynamics of the QCD matter produced in heavy ion 
collisions, for which parton energy loss is considered to be a sensitive probe. 

Prehminary RHIC data on single inclusive electrons in Au-Au collisions [15,16] 
extend over a wider transverse momentum range, up to ~ 10 GeV, than 
previously published results [17]. The purpose of this letter is to assess to what 
extent these spectra and their azimuthal asymmetry [18, 19] are in agreement 
with our current understanding of medium-induced parton energy loss, and 
how they can help us to access in more detail the microscopic mechanisms un- 
derlying high-pr hadron suppression. To this end, we determine the sensitivity 
of the nuclear modification factor for single electrons -R^a(^'t) on uncertainties 
in the b- and c-benchmark cross sections, we quantify the possible contamina- 
tion of Raa{Pt) from other hard processes (namely Drell-Yan), and we study 
its dependence on those properties of the nuclear matter, which are expected 
to induce the energy degradation of hard parent partons. 



2. Uncertainties of the proton— proton benchmark 

Heavy quark production in proton-proton collisions can be calculated in per- 
turbative QCD via the collinear factorization approach [20-22]. In Fig. 1, 
we show the resulting electron spectrum for y^J/vjv = 200 GeV pp collisions 
at RHIC. This spectrum was obtained [23] by complementing a fixed-order 
plus next-to-leading-log-resummed (FONLL [24,25]) heavy quark pt distribu- 
tion with proper fragmentation functions, describing the hadronization into 
heavy hadrons, and with the subsequent decay of the heavy hadrons into elec- 
trons. The Drell-Yan contribution, also shown in Fig. 1, is obtained using 
PYTHIA [26] , with a Drell-Yan cross section rescaled such that the e"*" e~ 
invariant mass distribution matches the one given by the NLO calculation 
in Ref. [27]. The curves in Fig. 1 reflect central values of a calculation, in 
which the Drell-Yan contribution to the single electron spectrum is about 
10% at py = 10 GeV. This implies that neglecting the Drell-Yan contribution 
in computations of the nuclear modification factor of single electrons, one re- 
duces R^A by up to 0.1. However, the theoretical uncertainties are large. At 
smaller transverse momenta p^, on which we focus in this paper, Drell-Yan 
becomes even smaller (see Fig. 1), and will be neglected in what follows. 
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Fig. 1. (Color online) Contributions from semi-leptonic heavy quark decays 
and from Drell-Yan pairs to the single inclusive electron spectrum in proton-proton 
collisions at -^s = 200 GeV, calculated with NLO accuracy. 



In Fig. 2a, we compare the perturbative FONLL result for the single electron 
spectrum in ^/s = 200 GeV pp collisions [23] to data from the PHENIX col- 
laboration at RHIC [28]. Upper and lower lines of the different decay con- 
tributions indicate theoretical uncertainties: they were estimated by vary- 
ing the factorization and renormalization scales independently in the range 
mT/2 < fiF,R < 2mT with the constraint 0.5 < fip/fJ'R < 2. In addition, 
we have varied quark masses around their central values = 1.5 GeV and 
rrih = 4.75 GeV over the range 1.3 < rric < 1.7 GeV, 4.5 < mt, < 5.0 GeV. 
These effects were then added in quadrature. The uncertainties related to the 
fragmentation picture, intrinsically ambiguous in the pt ^ ttl region, were not 
exphcitly considered here, as they are not expected to be very large (see for 
instance [29]). Within errors, the comparison between FONLL predictions and 
experiment is fair, see Fig. 2a. However, the central value of the FONLL cal- 
culation under-estimates the central value of the measured electron yields by 
a factor 2-3. Thus, an ad hoc increase of the charm cross section by a factor 
of three, large but still at the borderline of the above-mentioned uncertain- 
ties, would bring the central values of experiment and theory into agreement 
(data not shown). This illustrates that there is only limited control over the 
proton-proton baseline, on top of which one aims at establishing medium ef- 
fects. Another illustration of this is that the 6-quark decay contribution may 
start dominating over the c quark one at a transverse momentum as low as 
= 2.5 GeV or as high as p^ = 9.5 GeV (see Fig. 2). This translates into a 
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significant uncertainty in -R^^, as we discuss now. 
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Fig. 2. (Color online) (a) Left: Comparison of the FONLL calculation of single 
inclusive electrons [23] to data from pp collisions at ^/s = 200 GeV [28]. Upper 
and lower lines are estimates of theoretical uncertainties, obtained by varying scales 
and masses, for details see text, (b) Right: The nuclear modification factor i?^^ of 
electrons in central Au-Au collisions for an opacity of the produced QCD matter 
characterized by the time-averag ed BDMPS transport coefficient q = U GeV^/fm. 
The shaded band indicates the theoretical uncertainty of the perturbative baseline 
only. Red dashed and blue dotted curves show R%j^ for c-quark and 6-quark decay 
contributions, respectively. Data taken from Ref. [15,16]. 



3. The nuclear modification factor for single electrons in Au— Au 
collisions at RHIC 



The medium-induced suppression of high-p^ electrons in nucleus-nucleus col- 
lisions can be characterized by the nuclear modification factor R\j^{pt), which 
compares the production of electrons in AA collisions to an equivalent number 
of pp collisions, 
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Here, (N^-^) is the average number of inelastic nucleon-nucleon collisions in a 
given centrality class. We calculate R%y^ according to the model used success- 
fully for the description of suppressed high-p^ light-flavored hadron produc- 
tion [6,7]. To calculate d^N^^^'^/dpxdy, we start from the FONLL spectrum 
of final state heavy quarks in pp collisions. The nuclear modification of par- 
ton distribution functions is expected to be at most a 10% effect and can 
be safely neglected in our calculation [14]. We calculate the medium-induced 
energy degradation of the partonic spectrum by means of quenching weights 
for massive quarks [2,14], before fragmenting the energy-degraded partons in 
the vacuum according to FONLL fragmentation functions. This assumes that 
the interplay between quenching and hadronization can be neglected. Our 
treatment of medium effects is based on a realistic description of the collision 
geometry used to compute the in-medium path length and medium density 
on the parton-by-parton level [6,7]. 
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Fig. 3. (Color online) The nuclear modification factor for single electrons in 
central Au-Au collisions at RHIC. Curves indicate the suppression for different 
opacities of the produced matter. The shaded band indicates the theoretical uncer- 
tainty of the perturbative baseline for q = 14GeV^/fm. 

The effect of parton energy loss on single inclusive spectra depends on the color 
charge and mass of the parent parton, its in-medium path length L, and the 
time-averaged squared momentum transfer from the medium to the partonic 
projectile, which is characterized by the BDMPS transport coefficient q [30]. 
Previous studies have shown that for 4 < g < 14GeV^/fm, this model of par- 
ton energy loss accounts for the strength and approximate py-independence 
of light hadron suppression at RHIC [6,7]. In Fig. 2b, we show the nuclear 
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modification factor (3.1) for q ~ 14GeV^/fm. For this large value of q, a cal- 
culation anchored on the FONLL perturbative baseline tends to over-estimate 
R\a- However, due to the sizable theoretical and experimental errors, claims 
about an inconsistency between theory and experiment are not supported by 
Fig. 2b. We observe that an experimental separation of h- and c- decay elec- 
trons would strongly enhance the sensitivity to the mass hierarchy of parton 
energy loss. Indeed, since the mass dependence of charm quark energy loss 
vanishes as a function of (mc/p^)^, the suppression of electrons from charm 
decays is comparable to that of light-flavored hadrons. In contrast, electrons 
from 6-decays are much less suppressed due to the larger mass of their quark 
parents, see Fig. 2b and Ref. [31]. 

In Fig. 3, we compare the same data on the nuclear modification factor to 
central values of the parton energy loss model, calculated for different values 
of the transport coefficient q. By comparing Fig. 3 with the corresponding 
results for the nuclear modification factor i?^^ for fight hadrons [6,7], we find 
that for the experimentally favored range q = A — 14GcV^/fm, the difference 
between R\j^ and i?^^ is roughly 0.2. In general, the theoretical uncertainty 
in calculating the partonic baseline spectrum is comparable to the model- 
intrinsic uncertainty in determining q. For instance, the central value for q ~ 
10 GeV^/fm lies within the perturbative uncertainty band of ^ = 14 GeV^/fm. 
We also note that as q increases, the suppression gradually saturates. This 
is due to the fact that at high opacity, particle production becomes surface 
biased, thereby decreasing the sensitivity of the nuclear modification factor on 
q [7,32]. Since heavy quarks lose less energy than light quarks or gluons, this 
surface bias becomes important at higher values of q. Thus, if experimental and 
theoretical errors can be controlled, the measurement of the decay products 
of heavy quarks docs have the potential to improve the characterization of the 
very opaque medium observed at RHIC [31]. 

Galculating in 'minimum bias' (0-80 %) Au-Au colfisions the angular de- 
pendence of the yield of heavy quark decay electrons with respect to the 
reaction plane, we can determine the azimuthal asymmetry, expected from 
parton energy loss in a spatially asymmetric medium [33,34]. This asymmetry 
is characterized by the so-called elfiptic flow, which is the second harmonic 
coefficient in the Fourier decomposition of the angular dependence with 
respect to the orientation of the reaction plane, dN^^^*^ jpT dpr d(j) cx 
1 + 2v2{j)t) cos [2(0 — ^_r)] [35]. At low transverse momenta, other collective 
phenomena are expected to dominate for heavy quarks and their decay 
electrons [36-41]. However, for pf, > 2 GeV, the azimuthal dependence of the 
energy loss of heavy quarks may be the dominant contribution to the observed 
asymmetry. As seen in Fig. 4, preliminary data in this kinematic regime show 
significant uncertainties, and our calculation favors a non-vanishing but small 
value for the high-p^^ azimuthal asymmetry of electrons. In the presence of 
collective transverse flow effects, the BDMPS transport coefficient may be 
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sensitive to a combination of energy density and directed collective motion. 
This is expected to further increase high-p^ f| [34,42]. Thus, our present 
study does not allow us to exclude values of f 2 which are somewhat larger 
than V2 = 0.05, but it strongly disfavors V2 > 0.1. 
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Fig. 4. (Color online) The azimuthal asymmetry of the single inclusive electron 
spectrum in Au-Au collisions at 0-80 % centrality, compared to our model of parton 
energy loss. Data taken from Refs. [18,19]. 



4. Conclusions 



Our study indicates that for the parameter range g = 4 — 14GeV^/fm of 
the BDMPS transport coefficient favored by the suppression of light hadron 
spectra at RHIC, the nuclear modification factor R^a of single electrons is 
~ 0.2 larger than that of light-flavored hadrons. Thus, single electron spectra 
are in principle sensitive to the expected mass hierarchy of parton energy loss. 
Moreover, since massive quarks lose less energy than light ones in the medium, 
their spectrum starts being dominated by surface emission for higher medium 
densities only. Thus, single electron spectra may increase the experimental 
sensitivity to the properties of the produced matter in the region of high 
opacity. 

However, to exploit these opportunities of single electron measurements in 
nucleus-nucleus collisions demands a significantly better control of the vari- 
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ous sources of uncertainties. The preliminary data available so far neither al- 
low us to support claims of an inconsistency between theory and experiment, 
nor do they support yet the theoretically expected mass hierarchy of nuclear 
modification factors. Here, we have contributed to the generally needed as- 
sessment of uncertainties by characterizing the theoretical ones. In particular, 
we have established that uncertainties in the perturbative baseline are signif- 
icant, and should not be neglected in calculations of the nuclear modification 
factor. If future experiments succeeded in disentangling the b- and c-decay 
contributions (e.g. by reconstructing the displaced decay vertices), this would 
remove a significant source of theoretical uncertainty, dramatically enhancing 
the experimental sensitivity to the mass hierarchy of parton energy loss, as 
can be seen from Fig. 2b. As an aside, we note that single electron spectra are 
not the only possibility to characterize the mass dependence of parton energy 
loss. Decay muons may provide another, as yet unexplored possibility. Also, 
topologically reconstructed D- and 5-mesons provide an alternative and much 
more versatile access to the same information, which is particularly promising 
at the LHC [14,43]. 

We finally emphasize that, although preliminary RHIC data favor high values 
{q = 14GeV /fm) of the BDMPS transport coefficient, significantly smaller 
densities given by ^ = 10 GeV^/fm reproduce the data almost as satisfactorily, 
since their theoretical uncertainties largely overlap (see Fig. 3). The densities 
required in our calculation to obtain -Rj4^(p^ ~ lOGeV) = 0.4 are not larger 
than those used previously [6, 7] for the characterization of light hadron sup- 
pression. 
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